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to an external magnetic field. The material used is PDMS reinforced with nickel 
coated carbon fibres. Thanks to the nickel functionalisation, the fibre orientation 
can be achieved by embedding the viscous solution into a low external magnetic field 
(< 0.2 T). It is shown that both mechanical and magnetic properties can be controlled 
by tailoring the material anisotropy through properly orientating the reinforcing fibres 
in the pre-curing phase. The large strain behaviour is investigated by tensile testing 
up to 60% of deformation and shows a strong dependence on the fibre orientation. The 
magnetic properties are investigated by placing beam-like specimens into a uniform 
magnetic field. The results show a multistable behaviour with a transition from a 
bending-only deformed configuration for the 0° fibres specimen, to a twisting only 
configuration, achieved for fibres at 90° whereas all the intermediate angles show both 
bending and twisting. This behaviour is accurately captured by the large rotations 
beam model introduced. Such an actuator can be used in all applications which 
require fast response times and large strain. 
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I. INTRODUCTION 

To enhance the mechanical properties of elastomer materials, reinforcements are often 
used. The use of discontinuous fibres is often motivated by the increased surface area and 
therefore the increased surface interaction with the matrix, when compared to spherical 
inclusions (such as Carbon black). However they are usually randomly dispersed to give an 
overall isotropic behaviour, except in those cases when the manufacturing process imparts 
some anisotropy (e.g. extrusion). 

However, by controlling the orientation of the discontinuous fibres in an elastomeric ma¬ 
terial, is it possible to affect both the small and large strain behaviour of the material. 
Conventional composite materials often utilize this orientation with continuous fibres to 
counter many different problems (e.g. see forward swept divergent wing- 10 ), however addi¬ 
tional non-linear effects can be observed when using discontinuous reinforcements at large 
strain (>20%)P®. This provides additional design space as it is possible to control the rein¬ 
forcement at discrete locations in 3-dimensions, to affect the material performance at large 
and small strain. 

The challenges in controlled orientation of materials in solution have been studiecP®, 
and have included various manufacturing methods, for example, sheai®, eleetrieaf 14 and 
magnetic alignment 2 1 5 1 17 1 ; the latter of which has demonstrated controlled 3-dimensional 
reinforcement 5 . Significant research has been employed to understand the factors under¬ 
pinning the orientation of magnetic fibres in a solution, such as field strengtbP® field 
direction®, and particle morphologjJ®, in order to spatially control the reinforcements. This 
has delicate consequences upon the mechanical properties, but can also add a degree of 
multi-functionality. For instance, the aligned fibres can also provide control of other mate¬ 
rial properties such as the electrical conductivit y* 11 * 14 ! dielectric 1 ®, energy dissipation®, gas 
permeability 19 and thermal conductivity 3 . 

It is also possible to utilize the reinforcement and material properties to actuate a spec¬ 
imen of the material. Actuation can be achieved by a variety of different methods, such 
as pH sensitivity^®, swelling 6 ®, electrical stimulation 1 and magnetic attraction 12 . The key 
advantage of aligning magnetic reinforcements, and using that induced magnetic responsive¬ 
ness of the material, is that the actuation direction is controlled; both by the non-isotropic 
mechanical properties and the non-isotropic magnetic responsiveness. 
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In this study we take inspiration from nature, and apply some of these concepts to the 
large strain behaviour of elastomer materials reinforced with discontinuous fibreP, although 
the work is equally applicable to any material with a viscous pre-cure phase. Rectangular 
PDMS specimens with different fibre angles were produced and tested to evaluate the influ¬ 
ence of the fibre orientation upon mechanical and magnetic properties. It is shown that the 
specimens respond to a low magnetic field (< 0.2 T) with large shape variations. Instability 
effects were also observed and were accurately captured by the large rotations beam model 
introduced. 

Whilst the focus of this investigation has been on the use of a transversely isotropic 
configuration, an out-of-plane reinforcement 12 or variable angle configuration could be pro¬ 
duced by such a method. The use of a nickel coating increases the permittivity of the fibre 
and therefore allows the use of a relatively low magnetic field that can easily be generated 
through single or multiple neodymium magnets. The material produced shows the potential 
for not only improving the properties of elastomer materials in-plane, but also spatially; 
with consequences upon the mechanical behaviour and multi-functional capabilities. 


II. MATERIALS PREPARATION 

PDMS (Sylgard 184) is used without further modification. Nickel coated carbon fibres 
(40% Ni/C chopped to ~0.25mm, diameter 4.8 pm) are purchased from Marktek Inc. and 
used as received. NiC short fibre composites are made by direct mixing of the two com¬ 
pounds; this involves the addition of a carefully weighed amount of NiC fibres to a test tube 
of ~24g PDMS. This mixture is firstly agitated to partially homogenise the mixture, then 
sonicated for 1 hour to further separate agglomerated fibres, and then again agitated. The 
mixture is then left for 1 hour so that any remaining agglomerates settle to the bottom; this 
assures a homogeneous dispersion of fibres upon curing. The supernatant is then separated 
and added to curing dishes in lg quantities and curing agent is added at a ratio of 1:10, 
with the remaining PDMS residue and NiC sediment cured. 

Each curing dish is then positioned between four Neodymium N52 magnets in an approx¬ 
imate homogeneous magnetic field of ~0.075 Tesla (see Fig. [l]). The mixture is then allowed 
to cure at room temperature for 48 hours. The magnetic setup is assembled from plastic, 
and arranged so that an approximately homogeneous magnetic field is present on the curing 
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material. Preliminary tests are performed on the arrangement to ensure the homogeneity of 
the fibres post-cure, and that orientation of the fibres has occurred. The homogeneity of the 
magnetic field is of importance so that the fibres orient themselves parallel to the magnetic 
field due to the magnetic torque; if the magnetic field was non-homogeneous then the fibres 
would also translate and hence agglomerate close to the boundary of the petri dish. As such 
the resulting material would be non-homogeneous. 



Figure 1. (a) Magnetic setup for specimen curing. The specimen is suspended in the centre of 
four Neodymium N52 magnets in a homogeneous magnetic field, (b) The magnetic field intensity 
calculated in the FE environment FEMM 15 . 

Once cured, the specimens are cut from the petri dish at predefined fibre orientation 
angles between 0° and 90° degrees, 2 specimens cut from each petri dish with a dumbbell 
cutting tool (Fig. [2]). A total of 6 specimens are made at each angle 9 q =[0, 15, 30, 45, 
60, 75, 90]. A concentration of NiC of 6.0 wt% is calculated post-cure. This is calculated 
by burn-off of the PDMS mixture. The residue is added to a crucible and heated to 500°C 
for 6 hours, the remnants of this burn-off are then weighed and then compared to that of a 
pure NR burn-off. This gives the amount of NiC not used in the specimens; with the total 
amount of NiC known, the amount used in each specimen can be calculated. 

Prior to mechanical testing, an image of each specimen is captured by optical microscopy 
(see Fig. [2]). This micrograph of the specimen is post-processed so that a 60 mm x 40 mm 
section of the specimen is analysed to determine a representative set of data regarding the 
position, orientation and size of the particles in each specimen. Image software is used to 
create an image layer onto which the the particles are drawn on; this layer is then imported 
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Figure 2. Dumbbell test coupon and inset micrograph of a 60° specimen. Selected fibres highlighted 
in white. 

to the matlab Image Processing Toolbox. 

Figure [3^ & b show the typical distribution curves for a specimen. The orientation shows 
a significant preferential direction aligned to the direction of the magnetic field, although 
there is an underlying homogeneous dispersion of fibres that do not orient. Presumably 
the nickel coating is no longer present in these cases. There is also some distribution in 
the length of the specimen, which may have some contributing effects upon the ability to 
orientate in a magnetic field (see Reference^. 

III. EXPERIMENTAL RESULTS AND DISCUSSION 
A. Mechanical Properties 

All mechanical tests were performed at room temperature, with dumbbell specimens 
conforming to ASTM standard D1708, with a specimen thickness of approximately 0.5 mm. 
The nominal stress is determined as the ratio of the measured force to the original specimen 
cross-sectional area, whilst the strain is the ratio of the extension to the original distance and 
is calculated by digital image correlation (DIC). Uniaxial mechanical testing is performed 
using an Instron testing machine. The specimens are tested applying displacement control 
(10 mm/min) until rupture. 6 specimens at each angle were tested. 

The results of the large strain tensile tests are shown in Fig. [4] for 0°, 30°, 60° and 75°. 
The differences in the stress-strain curves highlight the different stress transfer mechanisms 
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Figure 3. Distributions extracted from a micrograph of a 0° Specimen, sampling >600 fibres. 
Average parameter values shown with standard deviations for 6 samples, (a) Fibre Orientation 
Distribution, p(Q) = exp(6 cos(2Q))/ (2 exp(b cos(29))d6) (b = 4.49 ±0.57). (b) Fibre Length 
Distribution, h(L) — ex p( ~( 2 ^)' > ) ( c = 0-068 ±0.011, d = 0.089 ±0.017). 

that take place during deformation and cause the large variation of the tangent modulus 
seen in the figure. When the fibres are oriented in the loading direction (0° and 30°) they are 
able to carry the majority of the load and a high modulus is achieved; on the contrary low 
moduli correspond to configurations where the fibres are almost orthogonal to the loading 
(60° and 75°). 

It also is interesting to observe the effects of the fibre rotation on this stress transfer 
mechanism, which can be evidenced when comparing the reinforcement at 0° and 30° in 
Fig. The initial modulus of 0° is higher, but at A = ~1.35 the rotation of the fibre 
towards the loading direction causes the 30° tangent modulus to decrease less and the stress 
to rise above the other. A similar occurrence is not observed between the 60° and 75° 
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Figure 4. (a) Large strain experimental results of four specimens compared to the transversely 
isotropic (TI) constitutive model derived from Eq. [l] (b) Stress-strain curves, showing confidence 
intervals of 90% derived from the stress-strain data of each angle. 

specimens, however this is most likely due to the relatively moderate change in modulus 
that a rotation of the fibres would cause at these angles (see Fig. [5]). 

These results can be effectively interpreted by comparison with the transversally isotropic 
hyperelastic model used in Reference. The main underlying hypothesis of the model was 
that the strain energy density of the reinforced material is separable into two contributions: 
one accounting for the overall stretch of the PDMS matrix and the other weighing the stretch 
along the stiff fibres. This model is based on the assumption that the strain energy stored in 
the fibres depends only on the fibres elongation 4 9 . The corresponding nominal stress (force 
over reference area) in the loading direction is given by 


S( A) = T 4 A 3 cos 2 (d 0 ) - T 2 A 3 + Ti A 3 - p A, 


(1) 


with 6*o being the initial orientation of the fibres, A := £/£q the stretch and dR, \l / 2 and d / 4 
constitutive functions defined as 


Ti — fj,, d/ 2 = a // d/ 4 = /31 — 


A 3 / 2 


( 2 ) 


(A 3 cos 2 (# 0 ) + sin 2 (# 0 )) 3/2 

The unknown pressure field p in Eq. ([!]) can be worked out by enforcing the plane stress 
condition (null stress in the through-thickness direction). Note that the constitutive param¬ 
eter (3 weighs the anisotropy of the specimen and in fact the isotropic Mooney-Rivlin model 
is recovered if f3 = 0. 
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Figure 5. Variation of the Young’s modulus in terms of fibre orientations fitted to Eq. [3] 
(/i = 1.579xl0 _1 MPa, a = 1.850, (3 = 2.878). Confidence intervals of 90% shown for each angle. 

For such a model, the expression of the effective tangent modulus (Young’s modulus) in 
the direction of the tensile loading is 

E(6 0 ) = — [16 + 16a + 5(3+ 8/3 cos (20 O ) + 3 (3 cos (40 o )l. (3) 

16 

Equation[3]is able to describe the dependence of the Young’s modulus on the fibre orientation 
well, as shown by the fitting in Fig. [5j E(0) has a minimum for fibres oriented at do = 65.9° 
as corroborated by the stress-strain data in Fig.|4j For this angle the fibre provides the least 
amount of resistance to the uniaxial forces imposed upon the specimen. On the contrary 
for 0° or 90°, the effective modulus has a local maxima (see Fig. [5]). The maximum at 0° 
being significantly larger than at 90°. However, the limiting value at 90° is interesting as 
it indicates the potential importance of the lateral constraint of the reinforcement in an 
incompressible system. 

At strains above ~30% the model begins to struggle to capture the highly non-linear 
behaviour observed in the experimental data in Fig. [4ji. Therefore a more complex form 
of the constitutive parameters Ti, + > and 'Fi may be considered but it would be out of 
the scope of this paper; in fact, as shown in the next section, the deformations obtained 
by magnetically actuating a thin specimen are rarely larger than 5% making this model 
adequate to describe the material behaviour. 

To further analyse the transversally isotropic nature of the produced composite materials, 
during tensile testing, microscopy is undertaken on identical specimens at strains values up 
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to ~55% (see Fig. |6j). This has allowed the rotation of individual fibres to be captured and 
compared to rotations predicted by the constitutively transversely isotropic model. By using 
Nanson’s formula (see, e.g., Ref. 16 ), it is possible to obtain the current orientation of the fibre 
9 in terms of the original angle 9q and the stretch A, i.e., 9 = — tan -1 [A -3 / 2 tan(0 o )]. The 
results of the fitting with this formula are shown in Fig.[7]for five specimens with initial angles 
0°, 15°, 45°, 75° and 90°. A remarkably close agreement between the rotation predicted by 
the model and the actual rotation of the fibres is apparent up to 55% of deformation, which 
suggest the effectiveness of the plane stress assumption applied in Eq. Q. 



Figure 6. Experimental micrographs of fibre rotation in the PDMS matrix at different stretch levels, 
for a specimen with initial angle of the fibres at 45°. 



Figure 7. Rotation of individual fibres compared to the predicted rotation of the constitutive model. 
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Figure 8. Schematic representation of the beam model, (a) The undeformed beam showing the 
reinforcement angle in the x-z plane, do- (b) The twisting angle, ip. (c) The bending angle, i). 

B. Magnetic Actuation 

The actuation capabilities of the material were tested for each fibre orientation by cutting 
a rectangular shaped specimen and suspending it between the plates of an electromagnet in a 
cantilever configuration (as sketched in Fig. [8j) . The resulting magnetic field is homogeneous 
and its intensity can be controlled by varying the current absorbed by the electromagnet. 
Further to this, the dynamic response of the actuation was viewed by generating a rotating 
magnetic field at frequencies between 0 and 3Hz through cylindrical neodymium magnets; 
demonstrating its capabilities for fast response times. 

The advantages of having the fibres in a plane orthogonal to the magnetic field would 
be the possibility of having different actuator responses based on the fibre angle do- To 
investigate this behaviour we introduce a large rotations beam model for which the total 
energy is the sum of an elastic and a magnetic part, i.e., 

Etot = ht f (T e + T m )ds (4) 

Jo 

where h and t are the beam width and thickness respectively, and T e and T )n given by 

^'( s ) 2 + YfJ ^'( s ) 2 > = ~x(n.B) 2 . (5) 

Here the elastic energy is expressed in terms of the bending, d, and twisting, ip, an¬ 
gles and accounts for large rotations/displacements; 'F m depends upon the relative angle 
between the magnetic field B = {0,H y ,0} (here assumed to be directed along the y- 
axis) and the fibre orientation in the deformed configuration n = (sin(d 0 ) sin(d) sin(</?) + 
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cos(#o) cos($), cos(^o) sin(d) — sin($o) cos('d) sin(<p), sin(#o) cos(<p)} (see Fig. [8j). In Q s is 
the abscissa along the specimen length (0 < s < L) and x = Xa Ff/(2/io) where Xa represents 
the magnetic anisotropic susceptibility, // 0 the vacuum permittivity and Uf is the volume 
fraction of the fibres 12 . El and GJ are the bending and twisting rigidities, respectively. 

The deformed configuration of the beam specimen is studied by looking at the minima 
of the total energy; the model Q being nonlinear, only numerical solutions can be derived. 
However, a great deal of insight into the specimen behaviour can still be gained by taking 
a fourth order expansion of the total energy with respect to the maximum bending and 
twisting angles; as a result, one can obtain a closed form approximate expression of the 
stable and unstable equilibria configurations of the beam. This approximate expression of 
the total energy in terms of the maximum bending i) m and twisting ip m angles achieved at 
the free end of the beam is given by 


Erotai = ^B 2 y ^ m Lx cos 2 (#o) - ^B 2 ti 3 m Ltp m x sin(do) cos(do) + ^B 2 tf 2 m L(p 2 m x sin 2 (d 0 ) 

11 2 

- -B 2, d 2 m Lx cos 2 (# 0 ) - —B 2 y 'd m Lipl l x sin(d 0 ) cos(d 0 ) + -H 2 d m L<p m xsin(6> 0 ) cos(d 0 ) 


+ ^B 2 y Lip A mX sin 2 (0 o ) - ^L^xsin 2 (0 o ) + 


+ 


GJip: 


2 L 3 


( 6 ) 


The stability of the equilibria can be studied by looking at the second derivatives of 
Eq. (JgJ, i.e., the Hessian matrix of the system, which gives the following value of the critical 
magnetic field By Crit ^ 



3 El 


2 L 4 X (EI/GJ sin 2 (# 0 ) + cos 2 (d 0 )) ' 


( 7 ) 


Up to the critical value, By Crit \ the system has only one equilibrium position that corresponds 
to the undeformed configuration; as soon as B y > By llt) the initial configuration becomes 
unstable (local maximum of the energy) and two symmetric minima appear. In this situation 
the reinforcing fibres attempt to rotate parallel to the magnetic field. The presence of two 
symmetric minima is due to the fact that positive or negative bending/twisting angles have 
the same weights in the energy. 

The value of By Cnt ^ can be experimentally determined by slowly increasing the intensity 
of the magnetic field up to the point at which a sudden jump in the deformation of the 
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0.122 



Figure 9. Maximum bending angle in terms of the intensity of the magnetic field for a specimen 
with fibres at 9q = 30°. When the magnetic field overcomes ~0.122 T a sudden increase in the 
bending angle is observed. 


beam is observed. This effect is shown in Fig. [9] for a specimen with fibres at 6q = 30° 
for which it is seen that when the magnetic field overcomes ~0.122 T a sudden increase in 
the bending angle happens. This is clearly visible from the camera samples in the figure. 


The values of B^ rit ^ predicted by the model through Eq. 0 are shown in Fig. [l(j| against 
the experimental data for specimens with fibres at 0q = {0°, 15°, 30°, 45°}. For the other 
specimens the point of instability is less clear, and a more smooth increase in the bending 
and twisting angles with the magnetic field intensity is observed. This is partly due to 
the difficulty in measuring the twisting angle, which becomes the dominant actuation mode 
for higher angles, which is confused further by the fibres not being completely aligned but 
having a spread in the orientation as shown in Fig. (3}i. The fitting in Fig. [TO] was achieved 
with a value of the anisotropic susceptibility Xa = 1-21 x 10 -3 ; it should be noted that this 
value is about twenty times larger than the one reported in the literature for neat carbon 
fibred and is likely due to the nickel functionalisation of the fibre which make them easier 
to orient in low magnetic fields. 

The minimisation of (|6j) allows the equilibrium positions for each fibre orientation to be 


worked out as shown in Fig. 11; this latter was obtained with values of E(6*o) and G(6*o) taken 
from the mechanical characterisation carried out by the authors. At 0° the specimen exhibits 
bending only, whereas at 90° only twisting is observed. Fibre angles in between show both 
bending and twisting. The results of the experiments are shown in the same figure with 
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Figure 10. Critical magnetic field intensity By Crit ^ at which the undeformed configuration becomes 
unstable, assuming the model Q. (Parameters: E(0 q) and G(0o) taken from the results of the 
mechanical testing, I = 7.29xl0 -14 (m 4 ), J = 2.51xl0 -13 (m 4 ), L = 25xl0 -3 (m), \a = 1.21xl0 -3 
and Uf = 6 %). The green dots represents the results of the experiment carried out by the authors 
for specimens with fibres at 9q — {0°, 15°, 30°, 45°}. 



Figure 11. Maximum twisting and bending angles achieved for an imposed orthogonal magnetic 
field at the occurrence of instability, assuming the model (Parameters: E(0o) and G(0o) taken 
from the results of the mechanical testing, I = 7.29 x 10 -14 (m 4 ), J = 2.51 x 10 -13 (m 4 ), L = 25 x 10 -3 
(m), Xa = 1.21 xlO -3 and Uf = 6 %. The green dots represents the results of the experiment carried 
out by the authors for specimens with fibres at 6o — {0°, 15°, 30°, 45°, 90°}. 

green dots representing the bending angles of specimens with 0 O = {0°, 15°, 30°, 45°, 90°}. 
The corresponding twisting angles are not shown due to the large errors in their estimate; 
the width of the beam being much smaller than its length. 

Although a certain degree of dispersion in the fibre angle is present (see the distribution 
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function in Fig.[3j), the actuation of the specimens in the magnetic field confirm the behaviour 
predicted by the model. Indeed, as seen in Fig. [12} a 0° specimen exhibited only bending, 
whilst only twisting is observed at 90°. At angles between, either twisting or bending is 
observed as the dominant actuation mechanism; although it is still possible to observe both 


(see Fig. 12 3 & c). 



Figure 12. Static Actuation of selected specimens in a homogeneous magnetic field, (a) 0°, (b) 30°, 
(c) 45°, (d) 90°. Specimens of length ~25mm, width 7mm and thickness 1mm. 

IV. CONCLUSIONS 

We have studied the behaviour of magnetically responsive nickel-coated carbon fibres 
embedded in an PDMS elastomer material and oriented under a magnetic field to produce 
a transversely isotropic material. 

In this way, the material has mechanical properties that are controlled by the strength and 
direction of the magnetic field during the cure of the material to produce bespoke material 
properties. The fibre orientation causes both longitudinal and transverse stiffening effects 
that are observed both at small and large strains. 

The large strain macroscopic behaviour of the composite is also shown to be described 
by the transversely isotropic constitutive equation up to around 30%. In particular, the 
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rotation of the reinforcing fibres is closely matched to the experimental results. Above 30% 
strain, the model struggles to capture the non-linear behaviour of the material; likely due 
to the non-ideal interface between the fibres and the matrix and the consequent debonding 
of the fibres. 

Thanks to the nickel functionalisation of the fibres, the elastomer material is also respon¬ 
sive to low magnetic fields (<0.2 T) and is shown to actuate due to the reinforcing fibres 
attempting to align themselves parallel to the magnetic field lines. A beam-like configu¬ 
ration with fibres aligned within the plane of the beam has been studied. The resulting 
actuator shows a multistable behaviour for which, up to a critical value of the magnetic 
field, the undeformed configuration is stable and then becomes unstable for larger magnetic 
field intensities. Both twisting and bending behaviour is independently observed at different 
initial orientations of the fibre and predicted by the simple model introduced. 

The result of using magnetically responsive reinforcing fibres is that both mechanical and 
magnetic responses can be tailored, in a complimentary manner. This could potentially be 
useful when considering the significant design requirements of a micro-swimmer^, requiring a 
remotely actuated flexible material that still has significant stiffness to propel itself through 
a viscous fluid. The proper control of the direction of the external magnetic field could be 
used to accurately control the shape of the actuator in a robust manner which is a peculiarity 
of the proposed configuration. 
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